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Abstract 

A search for the A(4140) state in B + — > J/ip(f)K + decays is performed 
with 0.37 fb _1 of pp collisions at \fs = 7 TeV collected by the LHCb ex- 
periment. No evidence for this state is found, in 2.4 a disagreement with 
a measurement by CDF. An upper limit on its production rate is set, 
B{B+ -> X{4U0)K + ) x £(X(4140) -> J/i}<p)/B(B + J/^<j)K + ) < 0.07 at 90% 
confidence level. 
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In this article, results are presented from the search for the narrow X(4140) resonance 
decaying to J/ip<f) using B+ J/iP(f)K + events^ (J/ip <p -»■ K+K~), in a data 

sample corresponding to an integrated luminosity of 0.37 fb _1 collected in pp collisions at 
the LHC at ^/s = 7 TeV using the LHCb detector. The CDF collaboration reported a 
3.8 a evidence for the X(4140) state (also referred to as Y (4140) in the literature) in these 
decays using pp data collected at the Tevatron {^fs = 1.96 TeV) [l]. A preliminary update 
of the CDF analysis with 6.0 fb _1 reported 115 ± 12 B + — > J/ip(f)K + events and 19 ± 6 
V(4140) candidates leading to a statistical significance of more than 5 a pi. The mass anc 
width were determined to be 4143.4+^ ±0.6 MeV and 15.3^°^ 2 - 5 MeV, respectively^ 
The relative branching ratio was measured to be B(B + — > X(A1A0)K + ) x B(X(41A0) — > 
J/ip(f))/B(B + J/ilxpK + ) = 0.149 ± 0.039 ± 0.024. 

Charmonium states at this mass are expected to have much larger widths because of 
open flavour decay channels |3| . Thus, their decay rate into the J/ip <p> mode, which is near 
the kinematic threshold, should be small and unobservable. Therefore, the observation by 
CDF has triggered wide interest among model builders of exotic hadronic states. It has 
been suggested that the V(4140) resonance could be a molecular state j4- 10 , a tetraquark 
state 



11,12 , a hybrid state [13j|14j or a rescattering effect [15 16 . The Belle experiment 
found no evidence for the V(4140) state in the 77 — > J/ip(p process, which disfavoured the 
molecular interpretation [171. The CDF data also suggested that there could be a second 
state at a mass of 4274.41^ ± 1.9 MeV with a width of 32.3li5;3 ± 7.6 MeV (2). In this 
case, the event yield was 22 ±8 with 3.1<r significance. This observation has also received 



attention in the literature 18, 19 



The LHCb detector 20 is a single-arm forward spectrometer covering the pseudo- 
rapidity range 2 < rj < 5, designed for the study of particles containing b or c quarks. 
The detector includes a high precision tracking system consisting of a silicon-strip vertex 
detector surrounding the pp interaction region, a large-area silicon-strip detector located 
upstream of a dipole magnet with a bending power of about 4Tm, and three stations 
of silicon-strip detectors and straw drift-tubes placed downstream. The combined track- 
ing system has a momentum resolution Ap/p that varies from 0.4% at 5GeV/c to 0.6% 
at lOOGeV/c, and an impact parameter (IP) resolution of 20 \xm for tracks with high 
transverse momentum. Charged hadrons are identified using two ring-imaging Cherenkov 
detectors. Photon, electron and hadron candidates are identified by a calorimeter system 
consisting of scintillating-pad and pre-shower detectors, an electromagnetic calorimeter 
(ECAL) and a hadronic calorimeter (HCAL). Muons are identified by a muon system 
(MUON) composed of alternating layers of iron and multiwire proportional chambers. 
The MUON, ECAL and HCAL provide the capability of first-level hardware triggering. 
The single and dimuon hardware triggers provide good efficiency for B + — > J/ip<pK + , 
J/ip — > events. Events passing the hardware trigger are read out and sent to an 

event filter farm for further processing. Here, a software based two-stage trigger reduces 



the rate from 1 MHz to about 3 kHz. The most efficient software triggers 21 for this 



analysis require a charged track with transverse momentum (px) of more than 1.7 GeV 



1 Charge-conjugate states are implied in this paper. 
2 Units in which c — 1 are used. 



1 



(pr > 1-0 GeV if identified as muon) and with an IP to any primary pp- interaction vertex 
(PV) larger than 100 fim. A dimuon trigger requiring Pt(a<) > 0.5 GeV, large dimuon 
mass, M(yU + /i~) > 2.7 GeV, and with no IP requirement complements the single track 
triggers. At final stage, we either require a J/ip — » candidate with p^ > 1.5 GeV or 

a muon-track pair with significant IP. 

In the subsequent offline analysis, J/0 — > fi + fi~ candidates are selected with the 
following criteria: Pt(/-0 > 0.9 GeV, \ 2 P er degree of freedom of the two muons forming a 
common vertex, x^ tx (/i + /i _ )/ndf < 9, and a mass window 3.04 < M(/i + /i _ ) < 3.14 GeV. 
We then find K + K~K + combinations consistent with originating from a common vertex 
with Xvtx(-^ + -^ ~-^ + )/ndf < 9. Every charged track with p^ > 0.25 GeV, missing all 
PVs by at least 3 standard deviations (xi P {K) > 9) and classified more likely to be a 
kaon than a pion according to the particle identification system, is considered a kaon 
candidate. A five-track J/tpK + K~K + vertex is formed (xl tx (J/ip K + K ~K + )/ndi < 9). 
This B + candidate is required to have px > 4.0 GeV and a decay time as measured with 
respect to the PV of at least 0.25 ps. When more than one PV is reconstructed, the one 
that gives the smallest IP significance for the B + candidate is chosen. The invariant mass 
of a fi + fi~K + K~K + combination is evaluated after the muon pair is constrained to the 
J/ip mass, and all final state particles are constrained to a common vertex. 

Further background suppression is provided by a likelihood ratio. In the case of 
uncorrelated input variables this provides the most efficient discrimination between signal 
and background. The overall likelihood is a product of probability density functions, V(xi) 
(PDFs), for the four sensitive variables (xj): smallest Xip{K) among the kaon candidates, 
Xvtx^/V K + K~K + )/n<M, the pointing of the B + candidate to the closest primary vertex, 
Xjp(-B), and the cosine of the largest opening angle between the J/ip and kaon candidates 
in the plane transverse to the beam. The latter peaks towards +1 for the signal as the 
B + meson has a high transverse momentum. Backgrounds combining particles from two 
different B mesons peak at —1. Backgrounds including other random combinations are 
uniformly distributed. The signal PDFs, V S i g (xi), are obtained from the Monte Carlo 
simulation (MC) of B + J/ipK + K~K + decays. The background PDFs, V hkg ( x^J, are 
obtained from the data candidates with J/ip K + K~ K + invariant mass between 5.6 and 
6.4 GeV (far-sideband). A logarithm of the ratio of the signal and background PDFs is 
formed: DLL sig / bkg = — 2 ^^ln(P sig (xj)/7 :, bkg(a ; i))- A requirement on the log-likelihood 
ratio, DLL sig / bkg < —1, has been chosen by maximizing iV si g/ ^N sig + N hkg , where iV si g is 
the expected B + — > J/ipK + K^ K + signal yield and the iVbkg is the background yield in 
the B + peak region (±2.5 a). The absolute normalization of iV S i g and iVbkg comes from a 
fit to the J/ip<f)K invariant mass distribution with DLL sig /bkg < 0, while their dependence 
on the DLL si g/bkg requirement comes from the signal simulation and the far-sideband, 
respectively. 

The J/ip(j)K invariant mass distribution, with a requirement that at least one K + K~ 
combination has an invariant mass within ±15 MeV of the <fi mass, is shown in Fig. [TJ A 
fit to a Gaussian and a quadratic function in the range 5.1 — 5.5 GeV results in 346 ± 20 
B + events with a mass resolution of 5.2 ± 0.3 MeV. Alternatively requiring the invariant 
mass M(J/ip K + K~ K + ) to be within ±2.5 standard deviations of the observed B + peak 
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Figure 1: Mass distribution for B + — > J/ip(f)K + candidates in the data after the ±15 MeV <p 
mass requirement. The fit of a Gaussian signal with a quadratic background (dashed line) is 
superimposed (solid red line). 

position, we fit the M(K + K~) mass distribution (two combinations per event) using a 
binned maximum likelihood fit with a P-wave relativistic Breit- Wigner representing the 
0(1020) and a two-body phase-space distribution to represent combinatorial background, 
both convolved with a Gaussian mass resolution. The resonance width is fixed to the 
PDG value (4.26 MeV) [22]. The M(K + K~) mass distribution is displayed in Fig. [2] 
with the fit results overlaid. The fitted parameters are the yield, the mass (1019.3 ± 
0.2 MeV), the background yield and the mass resolution (1.4 ± 0.3 MeV). Replacing the 
two-body phase-space function by a third-order polynomial does not change the results. 
In order to subtract a non-5 contribution, we fit the M(K + K~) distribution from the B 
mass near-sidebands (from 4 to 14 standard deviations on either side) leaving only the 
yield and the two-body phase-space background yield as free parameters. After scaling 
to the signal region, this leads to 14 ± 3 background events. The background subtracted 
B + -»■ J/^<pK + yield (N B+ ^j^^ K +) is 382 ± 22 events. 

To search for the X(4140) state, we select events within ±15 MeV of the mass. Ac- 
cording to the fit to the M(K + K~) distribution this requirement is 85% efficient. Figure[3] 
shows the mass difference M(J/if)<f>) — M(J/ip) distribution (no J/0 or mass constraints 
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Figure 2: Invariant M(K + K~) mass distribution selecting B + — > J/tp K + K~ K + events in 
the ±2.5(7 region around the B + mass peak. The dashed line shows the two-body phase-space 
contribution. The small blue dotted 4> peak on top of it illustrates the amount of the background 
4> mesons estimated from the fit to the B + mass near-sidebands. 

have been used). No narrow structure is observed near the threshold. We employ the 
fit model used by CDF |2| to quantify the compatibility of the two measurements. The 
data are fitted with a spin-zero relativistic Breit-Wigner shape together with a three-body 
phase-space function (J 7 } 3 ^), both convolved with the detector resolution. The efficiency 
dependence is extracted from simulation (Fig. [4j and applied as a shape correction to 
the three-body phase-space and the Breit-Wigner function. The mass and width of the 
X(4140) peak are fixed to the central values obtained by the CDF collaboration. The 
mass-difference resolution was determined from the B + — > X(A140)K + simulation to be 
1.5 ±0.1 MeV. A binned maximum likelihood fit of the signal and background yields is 
shown in Fig. [3|a). The region above 1400 MeV is excluded since it is more likely to 
contain non B + — > J/ip <pK + backgrounds. By excluding also the region below 1030 MeV, 
where the three-body phase-space and signal yields are very small (0.5% and 3.5% of the 
yields included in the fit, respectively), we make our results less vulnerable to possible 
small contributions from the other sources. The fit shown in Fig. [3]^a) gives a X(4140) 
yield of 6.9 ± 4.9 events. Fitting the second state at a mass of 4274.44 MeV and with 
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a width of 32.3 MeV [2] does not affect the A" (4140) yield. Reflections of K<p reso- 
nances 23 24] and possible broad J/ip(j> resonances can also contribute near and under 



the narrow A(4140) resonance. To explore the sensitivity of our results to the assumed 
background shape, we also fit the data in the 1020 — 1400 MeV range with a quadratic 
function multiplied by the efficiency-corrected three-body phase-space factor (J 7 ^) to 
impose the kinematic threshold. The preferred value of the A(4140) yield is 0.6 events 
with a positive error of 7.1 events. This fit is shown in Fig. |3|( 




Figure 3: Distribution of the mass difference M(J/ip(p) — M(J/ip) for the B + — > J/ipcj)K + in 
the B + (±2.5 cr) and <f> (±15 MeV) mass windows. Fit of V(4140) signal on top of a smooth 
background is superimposed (solid red line). The dashed blue (dotted blue) line on top illustrates 
the expected V(4140) (V(4274)) signal yield from the CDF measurement (2). The top and 
bottom plots differ by the background function (dashed black line) used in the fit: (a) an 
efficiency-corrected three-body phase-space {F l g ); (b) a quadratic function multiplied by the 
efficiency-corrected three-body phase-space factor (J" 2 bkg ). The fit ranges are 1030-1400 and 
1020-1400 MeV, respectively. 



A similar fit was performed to simulated B + — > A(4140)A + data to estimate 
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Figure 4: Efficiency dependence on M(J/ip<fi) — M(J/ip) as determined from the simulation 
(points with error bars). The efficiency is normalized with respect to the efficiency of the </> 
signal fit to the B + — > J/il)(f)K + events distributed according to the phase-space model. A cubic 
polynomial was fitted to the simulated data (superimposed). 



the efficiency for this channel. The efficiency ratio between this fit and the sig- 
nal fit to the B + — > J/ip(f)K + events distributed according to the phase-space model, 
e(B + -> X(4140)K+,X(4140) J/tp 0) / 't{B + J/ip^K+), was determined to be 
0.62 ± 0.04 and includes the efficiency of the <fi mass window requirement. Using our 
B + — > J/if)(pK + yield multiplied by this efficiency ratio and by the CDF value for 
B(B + ->> X(4U0)K + )/B(B+ J/t/j<j)K + ) leads to a prediction that we should have 
observed 35 ± 9 ± 6 events, where the first uncertainty is statistical from the CDF data 
and the second includes both the CDF and LHCb systematic uncertainties. Given the B + 
yield and relative efficiency, our sensitivity to the X(4140) signal is a factor of two better 
than that of the CDF. The central value of this estimate is shown as a dashed line in 
Fig. [3j Taking the statistical and systematic errors from both experiments into account, 
our results disagree with the CDF observation by 2.4 a (2.7 cr) when using ^ kg (J-J^) 
background shapes. 

Since no evidence for the X(4140) state is found, we set an upper limit on its produc- 
tion. Using a Bayesian approach, we integrate the fit likelihood determined as a function of 
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the X(4140) yield and find an upper limit on the number of signal events of 16 (13) at 90% 
confidence level (CL) for the two background shapes. Dividing the least stringent limit on 
the signal yield by the B + -> J/ip(f)K + yield and e(B + ->- X(AU0)K + )/e(B + ->• J/ip(f)K + ) 
gives a limit on B(B + X(4140)if + ) x B(X(AU0) J/ip(f))/B(B + J/ip(j)K + ). The 
systematic uncertainty on e(B + — > X (4140) K + )/e(B + — > J/ip(pK + ) is 6%. This uncer- 
tainty includes the statistical error from the simulation as well as the observed differences 
in track reconstruction efficiency between the simulation and data measured with the 
inclusive J/ip —> /i + n~ signal. Fit systematics related to the detector resolution and the 
uncertainty in the shape of the efficiency dependence on the J/i/j cf> mass were also studied 
and found to be small. We multiply our limit by 1.06 and obtain at 90% CL 

B(B + -j- X(4U0)K+) x B(X(AU0) J/ijjcj)) 

B(B+ J/^j<f)K+) < 

We also set an upper limit on the X(4274) state suggested by the CDF collaboration 
[2^. The fit with J-~^ ks background shape gives 3.4^3 4 events at this mass. The fit with the 
J- 2 g background shape gives zero signal events with a positive error of 10. Integration 
of the fit likelihoods gives < 24 and < 20 events at 90% CL, respectively. The relative 
efficiency at this mass is e(B + X (4274) K+, X (4274) J/^(f))/e(B + J/^j(f)K + ) = 
0.86 ± 0.10. The least stringent limit on the signal events yields an upper limit of 

BjB+ -» X(4274)iT+) x g(X(4274) -> J/^0) 
B(B+ -> J/^(j)K+) 

at 90% CL, which includes the systematic uncertainty. CDF did not provide a mea- 
surement of this ratio of branching fractions. Assuming the efficiency is similar for the 
X(4274) and X(4140) resonances, their X(4274) event yield corresponds to B(B + — > 
X(A27A)K + ) x £(X(4274) J/iP(f))/B(B + -> J/ij(j)K + ) = 0.17 ±0.06 (statistical uncer- 
tainty only). Scaling to our data, we should have observed 53 ± 19 X(4274) events, which 
is illustrated in Fig. [3j 

In summary, the most sensitive search for the narrow X (4140) —> J/ip <j) state just above 
the kinematic threshold in B + — > J/ijj(j)K + decays has been performed using 0.37 fb _1 of 
data collected with the LHCb detector. We do not confirm the existence of such a state. 
Our results disagree at the 2.4 a level with the CDF measurement. An upper limit on 
B(B + ->■ X(4U0)K+) x B(X(AU0) -> J/ip 0))/ B(B + ->■ J/^(f)K + ) of < 0.07 at 90% CL 
is set. 



Acknowledgments 

We express our gratitude to our colleagues in the CERN accelerator departments for 
the excellent performance of the LHC. We thank the technical and administrative staff at 
CERN and at the LHCb institutes, and acknowledge support from the National Agencies: 
CAPES, CNPq, FAPERJ and FINEP (Brazil); CERN; NSFC (China); CNRS/IN2P3 
(France); BMBF, DFG, HGF and MPG (Germany); SFI (Ireland); INFN (Italy); FOM 



7 



and NWO (The Netherlands); SCSR (Poland); ANCS (Romania); MinES of Russia and 
Rosatom (Russia); MICINN, XuntaGal and GENCAT (Spain); SNSF and SER (Switzer- 
land); NAS Ukraine (Ukraine); STFC (United Kingdom); NSF (USA). We also acknowl- 
edge the support received from the ERC under FP7 and the Region Auvergne. 



References 

[1] CDF collaboration, T. Aaltonen et ai, Evidence for a narrow near-threshold structure 
in the J/ip(f) mass spectrum in B + — > J/ip(f)K + decays, [Phys. Rev. Lett. 102 ( 2009) 
|242002| |arXiv : 0903 . 2229| 

[2] CDF collaboration, T. Aaltonen et ai, Observation of the y(4140) structure in the 
J/ipifi mass spectrum in B — > J/ipipK decays, (arXiv: 1101 . 6058, 



[3] N. Brambilla et at, Heavy quarkonium: progress, puzzles, and opportunities, Eur. 



|Phys. J. C71 (2011) 1534 


,|arXiv: 1010.5827. 


X. Liu and S.-L. Zhu, F(4143) is probably a 


Rev. D80 (2009) 017502] 


arXiv: 0903. 2529. 



[5] T. Branz, T. Gutsche, and V. E. Lyubovitskij, Hadronic molecule structure of the 
F(3940) and F(4140), |Phys. Rev. D80 (2009) 054019| |arXiv : 0903 75424] 

[6] R. M. Albuquerque, M. E. Bracco, and M. Nielsen, A QCD sum rule calculation for 



the F (4140) narrow structure, jPhys. Lett. B678 (2009) 186[ |arXiv : 090375 540 



[7] G.-J. Ding, Possible molecular states of D*D* system and F(4140), |Eur. Phys. J. 



[C64 (2009) 297} |arXiv:0904. 1782| 

J.-R. Zhang and M.-Q. Huang, (Qs)^(Qs)^ molecular states from QCD sum rules: 
a view on Y (4140), |J. Phys. G37 (2010) 025005] |arXiv:0905.4178| 

[9] X. Liu and H.-W. Ke, The line shape of the radiative open-charm decay of Y (4140) 
and F (3930), Phys. Rev. D80 (2009) 034009[ |arXiv : 0907 . 1349| 

[10] Z.-G. Wang, Z.-C. Liu, and X.-H. Zhang, Analysis of the Y(A1A0) and related molec- 
ular states with QCD sum rules, |Eur. Phys. J. C64 (2009) 373| |arXiv:0907. 1467| 



[11] F. Stancu, Can Y(4140) be a cess tetraquark ?, J. Phys. G37 (2010) 075 017 
larXiv: 0906. 24851 

[12] N. V. Drenska, R. Faccini, and A. D. Polosa, Exotic hadrons with hidden charm and 



strangeness, Phys. Rev. D79 (2009) 077502, arXiv: 0902. 2803 



[13] N. Mahajan, F(4140): possible options, Ph ys. Lett. B679 (2009) 228 
larXiv: 0903. 31071 



S 



[14] Z.-G. Wang, Analysis of the F(4140) with QCD sum rules, |Eur. Phys. J. C63~2 009) 
[TT5] [arXiv : 0903 .5200] 



[15] X. Liu, The hidden charm decay of F(4140) by the rescattering mechanism, Phys. 
|Lett. B680 (2009) L37l |arXiv: 0904.0 136[ 



[16] D. V. Bugg, Explanation for Y ■(4140) and X(3915), arXiv: 1103.5363 



[17] Belle collaboration, C. Shen et ai, Evidence for a new resonance and search for 
the F(4140) in the 77 -+ (pJ/tp process, |Phys. Rev. Lett. 104 (2010) 1120041 
larXiv: 0912. 23831 

[18] J. He and X. Liu, The line shape of the open-charm radiative and pionic decays of 



y(4274) as a molecular charmonium, arXiv : 1 102 . 1 127 



[19] S. I. Finazzo, M. Nielsen, and X. Liu, QCD sum rule calculation for the charmonium- 



like structures in the J/ip(j) and J/ijjuj invariant mass spectra, |Phys. Lett. B 701 



|(2011) 101[ |arXiv: 110272 347 



[20] LH Cb collaborat ion, A. A. Alves Jr. et al, The LHCb detector at the LHC, |JINST| 
|3 (2008) S08005] 

[21] LHCb collaboration, V. Gligorov, C. Thomas, and M. Williams, The HLT inclusive 
B triggers, |LHCb-PUB-2011-016| 



[22] Particle Data Group, K. Nakamura et al., Review of particle physics, J. Phys. G37 
|(2010) 07502T] 

[23] D. Frame et al., A spin parity analysis of the <pK + system produced in the reaction 
K + p <j)K + p, (f) K + K~ at 13-GeV/c, |Nucl. Phys. B276 (1986)~667] 

[24] Bari-Birmingham-CERN-Milan-Paris-Pavia collaboration, T. Armstrong et al., A 
partial wave analysis of the K~<fi system produced in the reaction K~p — > K + K~K~p 
at 18.5-GeV/c, \md. Phys. B221 (1983) Tj 



9 



